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Abstract. We explore inflationary trajectories within randomly-generated two-dimensional 
potentials, considered as a toy model of the string landscape. Both the background and 
perturbation equations are solved numerically, the latter using the two-field formalism of 
Peterson and Tegmark which fully incorporates the effect of isocurvature perturbations. Suf- 
ficient inflation is a rare event, occurring for only roughly one in 10 5 potentials. For models 
generating sufficient inflation, we find that the majority of runs satisfy current constraints 
from WMAP. The scalar spectral index is less than 1 in all runs. The tensor-to-scalar ratio 
is below the current limit, while typically large enough to be detected by next-generation 
CMB experiments and perhaps also by Planck. In many cases the inflationary consistency 
equation is broken by the effect of isocurvature modes. 
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1 Introduction 

Historically researchers have hoped that there is a fundamental physical theory uniquely 
describing what we observe. More than that, many have hoped that it has a signature 
allowing us to distinguish it from any other theory we might have imagined along the way. 
While string theory may or may not be the answer to all our questions, what is exciting is 
that for the first time it gives us a theory that might be (see ref. [1] for a nice review of string 
theory in cosmology). But string theory has surprised those who work with it time and time 
again. First there was one string theory, then there were five, then it turned out there was 
a continuum of theories [2, 3]! 

How are we to interpret this? Susskind coined the term landscape [3]. He described 
what he called the megaverse (now more commonly referred to as the multiverse), where 
in the low-energy approximation, different regions may be characterized by the values of a 
large number of scalar fields. The consequence of this is that we have some very complicated 
potential V(4>i, ...,(/>£)), with a very large number of minima, each corresponding to a different 
metastable vacuum energy. This implies that instead of trying to predict the values of 
observables, we should be trying to predict probability distributions. 

An early study of the possible consequences of this landscape picture for inflation was 
carried out by Tegmark [4], who generated a large number of random one-dimensional po- 
tentials and explored the inflationary outcomes. However a one-dimensional approach gives 
a very limited view as compared to the possible dynamics of the landscape, in particular 
limiting the effect of the choice of initial conditions and restricting the observable outcome to 
adiabatic perturbations. In this article we take further steps towards a more realistic rendi- 
tion of the landscape, while remaining quite rudimentary, by carrying out a similar analysis 
in a two-dimensional field space. This broadens the effect of initial conditions, as there are 
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now a family of possible trajectories passing through each point in field space, and permits 
isocurvature perturbations which can modify the form of the late-time adiabatic spectrum. 
Our aim is to characterize the spread of observational predictions for such models. In fu- 
ture work we will extend further to a D-dimensional landscape, which adds computational 
complexity but no new issues of principle, unlike the extension from one to two dimensions. 

2 Approach 

2.1 Where we intend to explore 

Ideally we would like to explore the landscape potential V(</>i, </>z>) directly but sadly 
the explicit form is currently unknown. So to be getting on with we construct an artificial 
potential. The hope is we can explore how certain characteristics expected to be manifest 
in the true potential give rise to a probability distribution for observables. Adapting the 
approach taken in ref. [4], we define a potential 

V(4>x,..,4> D ) = m\f (—,.., — — ] , 

where / is a well-behaved dimensionless function and m v and m^, (i G \..D), are the 
characteristic vertical and horizontal mass scales. Ultimately m v is to be adjusted to give 
the correct amplitude of observed perturbations. Adjusting this mass does not affect anything 
other than the distribution of vacuum energies, so for now we will not concern ourselves with 
it. 

In ref. [4] the case where D = 1 was extensively investigated. We wish to move to higher 
D as it allows for a broader range of behaviour, particularly isocurvature perturbations which 
will be one of the main focuses of this paper. In this paper we only take the step of increasing 
D from 1 to 2; however we believe this is the biggest step as there are no further qualitative 
differences to go beyond D = 2. We reserve analysis of that case to future work. 1 Table 1 
summarizes some of the changes as the number of scalar fields is increased. 

To construct our landscape, following an approach similar to ref. [4] we use a random 
function of the form 

m n 
j=l k=l 

where in practice we truncate the series at the values m = n = 5 and the Fourier coefficients 
a j,k, bj t k are independent Gaussian random variables with zero mean and standard deviation 

a = e ~0' 2 +fc 2 )/2" ; (2.3) 

With this form, the potentials we simulate are periodic with periodicity scale 27rmh p and 
we can only expect reasonable results if the evolution spans a distance in the x—y plane less 
than the periodicity of the function. For our choice of the horizontal mass m ni , which we are 
about to discuss, the periodicity is large enough to have no effect. 

The precise form of our potential clearly has no theoretical motivation and as such we 
are free to tinker with it as we please. We are interested in how features of the landscape affect 

An exception may be the distribution of the vacua, which changes in an interesting way both classically 
and particularly with regard to tunnelling [5-7] . 
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Issues 


Requires very flat 
potentials that may 
be hard to realize in 
fundamental theory 


Initial Conditions 


Initial conditions 



Table 1. Summary of some of the differences between single scalar field models and more scalar 
fields. 



the evolution of the power spectrum and observables, so what we need to consider is on what 
scale these features occur. If any, it is this quality of our landscape which should be motivated 
by fundamental theory. Our choice of potential allows us to adjust this in three ways; the 
truncation number, the standard deviation of the Fourier coefficients, and the choice of m^. . 
The truncation number and tuning of the standard deviation are specific to our potential, 
but the mass scale is a rather more general feature of potentials. Thus, to try and focus on 
this sense of scale, we fix the deviation and truncation number as stated above and discuss 
our options in terms of adjusting the horizontal mass with respect to some reference mass, 
taken to be the reduced Planck mass Mp\. If we think of features in the landscape as being 
anything that can cause a change in direction of the inflationary trajectory, then we are 
essentially asking how many of these features we expect the trajectory to encounter during 
an evolution giving rise to a sufficiently large number of e-folds of inflation. There are many 
interesting mass scales that we will not be attempting to explore here, but crudely speaking 
they fall into four categories: 

m-h; 3> Mp\: This mass choice is poorly motivated by theory but was nonetheless explored 
as one of many cases for single scalar field in ref. [4]. More generally speaking, this 
case corresponds to a trajectory on a nearly flat, almost featureless potential. The 
advantage of this is that it is highly predictive since all trajectories look pretty much 
the same and in the single-field case can readily be in agreement with observation. 
Also, it is easy to get lots of inflation. 

m-h; -C Mp\: This mass choice corresponds to something akin to an egg box. 2 The main issue 
with this choice is that it becomes difficult to achieve sufficient inflation without getting 
eternal inflation, since the trajectory will almost always roll straight into a minimum. 
That said, one might imagine that if the number of scalar fields was sufficiently large 
then the chance of getting sufficient inflation would increase. This is a particularly 
important effect if tunnelling is taken into consideration [5-7], but to keep things simple 
we will not be doing so in this paper. 

2 Allowing a high truncation number when rrih i > Mpi introduces small-scale power into the potential to 
similar effect. 
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m-h; -C Mpi and 77%. 3> Mpp While our choice of potential would require additional terms 
to investigate this scenario, we are referring to the sort of case where there is more 
than one kinematically significant scale. An example of this kind of situation was 
investigated in refs. [8, 9]. The setup they considered could be imagined as a sort of 
multi-dimensional version of a board with nails stuck in it. The trajectory has a slow- 
roll drift velocity with what the authors describe as a brownian motion imposed on 
top of this. Depending on the scales involved, this could lead to interesting features in 
the power spectra and, for a given suitable background evolution, the extra distance 
covered due to the random motion will increase the number of e-folds. 

rrih; ~ Mp\: Finally we have the beginners' ski slope, an example of which is shown in Fig. 1. 
This is what we will be investigating. There are no jumps and the features are gentle 
so an advanced skier or snowboarder would probably be rather bored in our landscape, 
but for inflation we feel this is an interesting scale on which to begin our exploration. 
This mass scale is well motivated by theory and it also gives rise to quite a broad range 
of behaviour, since in order to obtain a sufficient number of e-folds of inflation, the 
trajectory generally has to take a non-straight path. The downside is that the broad 
range of behaviour will make the model less predictive but then again, it is interesting to 
see how robust the values of certain parameters are under such variations. Also, while 
the variability in a two scalar field model may be large, one can easily imagine that 
the deviation might decrease as more scalar fields are introduced. This tendency was 
seen in ref. [10] for Nflation models with random initial conditions on many uncoupled 
fields. 

2.2 How we explore it 

As we discussed, any model of inflation where the potential has multiple minima predicts a 
probability distribution for the cosmological parameters. We wish to compute this distribu- 
tion for the potential described above. To do this we perform essentially the same experiment 
that was performed in ref. [4]: 

1. Generate a random potential V(4>) and start at (f) = (0,0). 

2. If V(0, 0) < then reject model, otherwise evolve to find the field trajectory. 

3. If model gets stuck in eternal inflation, reject. 

4. Once the model stops inflating, if the number of e-folds of inflation N < 60 we reject 
as insufficient inflation occurred. Otherwise find the local minimum to calculate /9 V ao 
and if p va _ c < 0, reject. If N > 60 and p vac > 0, calculate observables. 

5. Repeat many times to obtain a statistical sample. 

6. (Change some assumptions and do it all again.) 

We start the evolution at V(0, 0) for practical purposes. Instead of viewing this as 
starting at the same position in many different potentials, since the statistics of our landscape 
are invariant under translation, one could equally view this as starting at random positions 
in one infinite potential. An alternative choice of initial conditions would be to start only at 
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Figure 1. An example of our landscape, shown on two scales together with the trajectory taken 
during inflation. The first scale corresponds to the periodicity of the function, and the second shows 
the features relevant to the inflationary trajectory. The "water level" indicates the height of the 
minima determining the vacuum energy. We see this is very close to V = 0; this, together with the 
fact that this trajectory gives a sufficient number of e- folds of inflation (in this case N = 173), means 
that we have an anthropically suitable bubble universe. 

maxima. This relates to the choice of measure, a problem we will discuss shortly but which 
in reality is an issue well beyond the scope of this paper. Step 2 is there since it would make 
no sense to do otherwise; generally we are not interested in the statistics of initial conditions 
that don't give rise to regions of space. Step 3 rejects trajectories that get stuck in minima 
that lead to eternal inflation. We will not be considering tunnelling in this paper, so once 
again we give no statistical weighting. 

Step 4 describes what we do when we finally do encounter a potentially viable trajectory. 
The rejections that take place here are more debatable since they are essentially anthropic 
rejections. In ref. [11] it was shown that anthropic arguments place an incredibly strong 
bound on the vacuum energy. This was and still is a very important result and the realization 
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that it sits so naturally with string theory is more than a little exciting. We do not explore 
this idea in our current work; indeed we ignore it since the extremely narrow permitted 
range for the present vacuum energy is not computationally accessible. Instead we place 
only an approximate lower bound, rejecting universes with a negative vacuum energy on 
the grounds that they will recollapse shortly after inflation ends. Technically very slightly 
negative vacuum energies are anthropically acceptable but this does not affect what we are 
interested in here. Rejecting inflation of less than 60 e- folds comes from arguing that this 
is a requirement for the formation of galaxies on a smooth cosmological background [4, 12]. 
Steps 5 and 6 follow trivially. 

In constructing our experiment in this way we have inadvertently chosen a measure. 
We do not wish to go into it in any depth here but it is important to realize that what is 
known as The Measure Problem [4, 13] directly affects this work. If we subscribe to Bayesian 
statistics then one can describe the current situation as follows: Define P(U) to be the 
probability of finding yourself in a particular region of space-time where inflation has ended 
which we now define as a bubble universe U, from now on simply referred to as a universe. 
Then define P(0) to be the probability of making an observation O in any bubble universe. 
What one defines as an observation changes the results dramatically. One could simply 
decide observation meant finding protons for instance. Here, observation corresponds to an 
observer measuring a particular set of cosmological parameters. The conditional probability 
of making an observation O in a universe U is P(0\U), called the likelihood. Meanwhile, 
the probability that we are in a universe U given that we have the observation O is P(U\0), 
called the posterior probability. Bayes' theorem then gives us 



Our experiment essentially makes an attempt at calculating P{0) but this can be written as 



where M is the multiverse, and dU is the measure. So as you can see we have issues. The 
probability of making an observation O given U is where anthropic arguments come in. For 
example, the probability of making an observation in a universe that can't form galaxies 
is, we guess, not very high. We also have the problem of counting infinities and this is the 
measure problem. We have possibly an infinite number of universes, infinite in extent and 
yet we need to count the probability of making observations in them. We do not wish to 
discuss this in any detail but we need to remember that our experiment, for the time being, 
uses a measure where initial conditions are weighted equally (step 1) to start with but then 
we give zero weighting to universes that are invalid on anthropic grounds (step 4). 

3 Calculating observables 

In this section we lay down the theoretical framework involved in calculating observables. 
Much of this section follows ref. [14], which gives an extremely thorough and comprehensive 
review of how to make two-field model calculations. We therefore only give a summary of 
the calculations involved, operating in a simplified setting where we assume the fields are 
canonically normalized and working to lowest order in the perturbations (since the effect of 



P{U\0) 



P{0\U)P{U) 

pJo) 



(2.4) 




(2.5) 
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higher-order terms would be much smaller than the uncertainty of our overall implementa- 
tion). First we introduce the unperturbed equations. In the context of our work these are the 
equations used to perform the experiment described in the previous section. Next we move 
on to the background kinematics. These will prove invaluable in our attempts to understand 
what aspects of the potential are important in affecting observables and how these character- 
istics extend to a larger number of scalar fields. We then move on to the perturbed equations 
which finally lead us to calculating observables. We use reduced Planck units (Mpi = 1) 
throughout and the vector notation <p = (cpi, .., <j)r>). In practice we will often simply have 
<j) = (0i,02) but where it will be helpful in later discussion, we write down explicitly how 
this extends to D scalar fields. 

3.1 Background equations 

We take the background spacetime to be flat, (3 + l)-dimensional, homogeneous and isotropic 
and thus described by the Robertson-Walker metric, 

ds 2 = -dt 2 + a 2 {t) [dx 2 + dy 2 + dz 2 ] , (3.1) 

where a(t) is the scale factor. We are investigating inflation driven by D scalar fields with 
canonical kinetic terms. The non-gravitational part of the action is 



r r i D 

S = J --Y^d^id^-Vi^.-Ao) 



gd A x, (3.2) 



i=l 

and thus we have the standard field equation 

4> + 3H<j> + vW = 0, (3.3) 

where H is the Hubble parameter, whose dynamics are found from the (0, 0) component of 
Einstein's equations, 

3H 2 = p. (3.4) 

From eq. (3.2) we also obtain 

P=\\^? + V , P= 1 -\<P\ 2 -V (3.5) 

3H 2 = ^\<j)\ 2 + V. (3.6) 

Working in terms of comoving time t is not particularly convenient for our current work. 
Instead it is much more helpful to make the transformation 

dN = Hdt, (3.7) 

where N is the logarithmic growth of the scale factor and represents the number of e-foldings 
of inflation. Using ./V is useful since it is directly related to observables; it also simplifies 
all the equations and renders them dimensionless, thus making it simpler to compare the 
relative size of terms. From now on we will use the notation 

s In < 3 ' 8) 



and so eq. (3.4) becomes 
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to represent differentiation with respect to N. 

To further simplify the equations we introduce the slow-roll parameter 



Inflation ends when a < <^=^ e > 1. If we combine the field equation (3.3) with 
the comoving time derivative of eq. (3.6) we see that the slow-roll parameter can also be 
interpreted in terms of the field speed; 



v = \4>'\ = V2e. (3.10) 

Making the transformation given in eq. (3.7) and substituting eq. (3.9), the field and Friedman 
equations simplify to 

6' + VhnV = ^— (3.11) 

3 — e 



and 



i V 

H = (3.12) 



respectively. 

Equation (3.11) together with eq. (3.9) give the evolution of the fields. To solve 
them we need to provide initial conditions. We choose these to be (ft = (0, 0) and eft' = 



3.2 Relating kinematics to the potential via the slow-roll slow-turn approxima- 
tions 

This subsection is predominantly a summary of ref. [14], although our case is simpler in that 
we are only concerned with canonical kinetic terms. We also make the extension to a larger 
number of scalar fields a little more explicit. 

We wish eventually to investigate the evolution of perturbations. As summarized in 
Table 1, when there is only one scalar field, the only type of perturbation that exists is density 
perturbations, generated by perturbing the field along the trajectory. When there are more 
scalar fields however, for every extra scalar field there is another direction perpendicular 
to the trajectory in which perturbations can also arise. The perturbations perpendicular 
to the trajectory correspond to isocurvature perturbations but these in turn can fuel the 
evolution of the density perturbations. Decomposing the perturbations in this way leads us 
to introduce a new basis. If we think of the old basis as B = {e±, ..,e£>} where e\ points 
in the direction </>i, e2 points in the direction <ft>2 and so on, then we can think of the new 
basis K = {eii,e_L 1 , -,e± n _ 1 } as a rotation of the old basis such that the first basis vector 
now points along the trajectory and we label it en. The other basis vectors will then point 
perpendicular to the trajectory and we label them e± i . We refer to this as the kinematic 
basis 3 and we denote the components of a general vector A and matrix M as, for instance 

A||=e r A , M u .=e\Me u . (3.13) 

3 For a fully kinematic basis, one could also define one of the perpendicular fields such that the direction 
of the turning of the field is entirely along it, but we choose not to do that here since our basis will make it 
clearer how to extend our statistics to D fields. 
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Using the kinematic basis we can now introduce the slow-roll slow-turn (SRST) ap- 
proximations. There is more than one way to interpret the well-known single-field slow-roll 
approximations and as discussed in ref. [14], this affects how the slow- roll approximations 
generalize to multiple scalar fields. With regard to how slow-roll extends, probably the 
minimal requirement is to say: 

• Expansion is nearly exponential. 

• Deviation from this expansion changes slowly. 

If we take this as our definition then our slow-roll conditions become 

1 



-y < i 



<C 1. 



(3.14) 



Although these are sufficient to guarantee the above requirements, it is helpful to have a 
slow-turn approximation which we define as 



< 1. 



(3.15) 



There are D — 1 of these slow-turn conditions and violation of any one of these conditions 
would render the trajectory no longer slowly turning, so one might prefer to instead write 



1 



D - 1 



«1, 



(3.16) 



but for our purpose we find the former more useful. If the field is both slowly rolling and 
slowly turning we say it satisfies SRST. 

Returning to eq. (3.11) we see that the left-hand side represents deviations from the 
SRST limit and so under SRST conditions we can write 



~ - V f In V 



and the Friedman equation becomes 



3H 2 



V. 



(3.17) 



(3.18) 



When trying to understand how various features of the potential affect observables, 
sometimes it is nice to interpret things in terms of the kinematics of the scalar field and at 
other times perhaps it is better to look directly at the underlying geometry. For this reason, 
amongst others, it is helpful to be able to approximately jump from one approach to the 
other. Direct from eq. (3.17) we have 



-|V f lnV\ 



and differentiating eq. (3.17) one obtains 

<j>" ~ -MV f In V, 

where M is the Hessian of In V, otherwise known as the mass matrix, 

M = V f VlnV. 



(3.19) 



(3.20) 



(3.21) 
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We therefore arrive at approximations relating the kinematic quantities directly to the po- 
tential: 



-Mi 



llll 



(3.22) 



termed the speed-up rate, and 



(3.23) 



known as the i-ih component turn rate. 

We don't look at the second-order equations here but we refer the reader to ref. [14] 
should they need them. 

3.3 The perturbed equations 

Continuing to adopt the same approach as that chosen in ref. [14] it can be shown that 
by using a multifield version of the Mukhanov-Sasaki variable [15], the field perturbations 
decouple from perturbations in the metric and thus we are able to focus solely on the curvature 
and isocurvature perturbations. The evolution of the field perturbations can be found by 
perturbing the equation of motion for the background. The standard result in Fourier space 
is 



+ mdcf) + ( k 

a 



H 

IP 



1 

H 



1 

H 



(3.24) 



where k is the comoving wavenumber. Substituting Eqs. (3.7), (3.9) and (3.11) one eventually 
reaches 



1 



b" + 56' + 



a 2 V 



M + 



(3 



(3.25) 



and we see that mode evolution is primarily governed by the mass matrix with small correc- 
tions. 

Rotating to the kinematic basis, we can separate eq. (3.25) to find evolution equations for 
the adiabatic (56^) and entropy modes (S(j)±). Starting with the adiabatic mode, projecting in 
the ey direction, then solving the resulting equation in the super-horizon limit (k/aH) 2 <C 1, 
one finds that the growing super-horizon adiabatic modes are described by 



'II 



D-l 



(3.26) 



or equivalently 



or when SRST holds 



D-l 



D-l 



Mil 



2 E M n 



(3.27) 



(3.28) 



So we see that the evolution of adiabatic modes can be inferred directly from the background 
kinematics. The first term of eq. (3.26) tells us that the faster the speed-up rate, the faster 



- 10 - 



the evolution of the mode. The second term tells us about the sourcing of the adiabatic 
modes from the entropy modes. We see that for a given size of entropy mode, the faster 
the turn rate, the more the adiabatic mode will be sourced by that entropy mode. For our 
experiment we will see there is a distribution associated with the likelihood of encountering a 
given turn rate. This means that as we extend to a larger number of scalar fields, it becomes 
increasingly likely that the fuelling from entropy modes will be a significant effect. 
Following the same procedure for the entropy modes we find that 

o — - + <tyj, = -fJ>l t • (3-29) 



where Hi is the effective entropy mass 



9 - e 



= M UU + ? — ^ ( ) (3.30) 



and so is well approximated by /j,±. ~ M± i ±.. The great news here is that eq. (3.29) tells us 
that each of the entropy modes evolves independently of the others. This means we can find 
their amplitude and determine the evolution of adiabatic modes without solving a massive 
set of fully coupled equations. We also see that the evolution of a given entropy mode is 
dominated by the curvature of the log of the potential along that entropic direction. This 
makes sense intuitively as when the curvature is positive, a given trajectory is stable in the 
sense that a perturbed trajectory will be redirected back onto the background trajectory. 
Conversely, when the curvature is negative, the trajectory is dispersive and perturbations 
will evolve, fuelling the growth of that entropy mode. 

From here we are in a position to discuss curvature and isocurvature perturbations. 
The curvature perturbation represents the perturbation in the curvature of constant-time 
hypersurfaces. In the comoving gauge we have the gauge-invariant quantity, R, which during 
inflation can be shown to be 



R=— -. (3.31) 

v 

Isocurvature perturbations represent relative fluctuations in the different fields that leave the 
total curvature unchanged and hence are related to entropy perturbations, here defined as 4 

S^ 6 ^. (3.32) 
v 

Using eq. (3.27) then gives us 

R = +2 Yl lv) Si ( 3 - 33 ) 

i ^ ' 

and so we see that the super-horizon evolution of the density perturbation is independent 
of the speed-up rate; instead it depends on the sum of the isocurvature perturbations, each 
multiplied by their corresponding turn rate. We also see that when there is no turn rate (in 
any of the directions), the single- field result that the density perturbation is conserved on 
super-horizon scales is recovered. 



4 This definition, following refs. [14, 16], is chosen so that the isocurvature mode typically has the same 
amplitude as the adiabatic one at horizon crossing, and differs from the alternative convention where the 
isocurvature perturbations are normalized like massless field perturbations. 
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M|m > 


M|m < 


M_l_l > 


Playoff 


Damping 


M ±± < 


Fueling 


Playoff 



Table 2. Origins of evolution of isocurvature modes. 



For the isocurvature perturbations we use the approach of ref. [16] and parameterize 
the isocurvature modes as 

Si = faSi. (3.34) 



We then find that (3 can be well approximated in the SRST limit as 

A ^ Mm, - M uu 



(3.35) 



So interpreting this from the geometrical point of view we see that the evolution of a given 
set of isocurvature modes depends on the integral of the difference between the curvature 
along the adiabatic and entropic directions in a very intuitive manner. As summarized in 
Table 2, negative curvature in the adiabatic direction and positive curvature in the entropic 
direction corresponds to maximum damping of isocurvature modes, while positive curvature 
in the adiabatic direction and negative curvature in the entropic direction corresponds to 
maximum fuelling. 

Plugging eq. (3.34) into eq. (3.33) we have the result 



D-l rN 



R ~ 



v 



(3.36) 



and the total isocurvature perturbation S is given by 

D-i D-l 



s = E s ^ E S^Wn-"^)*", 



(3.37) 



where * denotes that a given quantity is to be evaluated at horizon crossing. Continuing 
with the approach of ref. [16] it will be helpful to rewrite eqs. (3.36) and (3.37) as 



/ R \ (lT RSl -- - Trs D -i\ I R * \ 

Si 



T SSl 

\So-lJ \0 Tssn-J 
where the transfer functions Tr^ and Tss t are given by 



Su 
\Sd~i*J 



T RSl (N*,N) 



2^T ss .(N*,N)dN. 

N, V 



and 



(3.38) 



(3.39) 
(3.40) 



T SSi (N*,N) = e^ M H- M ^ dN 

respectively. To reduce computational effort we substitute the approximation of eq. (3.23) 
when calculating Trs (3.39). In the cases we tested we found this approximation to be more 
than satisfactory. 
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3.4 Power spectra and cosmological parameters 

We are now finally ready to discuss the evolution of the power spectra. For a general quantity 
X, the power spectra and cross spectra are defined as 

P x 5 3 (k-k)^^(X(k),X\k)} (3.41) 

and 

c xy s 3 (k-k) = ^<*(k),3;t(k)> (3.42) 

respectively. We do not go into details here but quantizing and solving the perturbed equa- 
tions (see e.g. ref. [14]) leads us to 

PR * " (fr*) 2 i [1 + 2{C " 1)6 " 2CM|II|] * ' (3 ' 43) 
H*\ 2 1 



HA 2 1 



and 

where C = 2 — In 2 — 7 ~ 0.7296 and 7 is the Euler-Mascheroni constant. For future use we 
also introduce the tensor power spectrum which takes the usual form 



P ^- 8 (^r) [1 + 2(C-1)4, (3.46) 

which is conserved for super- horizon modes. We used the above second-order expressions in 
our computation but as one would expect, there would have been almost no difference had 
we used first-order approximations. 

Applying the transfer matrix eq. (3.38) we finally arrive at one of our most important 
destinations, namely the power spectra at the end of inflation 

D-l 

Pr = Pr* + Y. ( 2T RS>C RS „ + T 2 RSt P s „) , (3.47) 

% 

D-l D-l 

Ps=Y, P s< = J2 T ss,Ps„ (3-48) 



and 

CRSt = Tss,Crs„ + TRSiTss.Ps,,- (3.49) 

We see most of our previous comments on super-horizon evolution apply in much the same 
way to the power spectra. It is now also easy to see one of the reasons why we chose to break 
our exploration down into two stages. The interesting super-horizon behaviour predominantly 
comes from the transfer functions. We have seen that the super-horizon modes associated 
with entropy perturbations evolve independently of one another, each separately affecting the 
evolution of the adiabatic power spectrum. This means that by making a statistical analysis 
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of one set of entropy perturbations, we can easily see how the statistics will generalize to a 
larger number of scalar fields. 

Now that we finally have the power spectra the hard work is over. From here it is 
straightforward to find approximations for what we consider the key observables, namely the 
tensor-to-scalar ratio, spectral indices, and the running. We will evaluate these 55 e-foldings 
before the end of inflation, taken to be the time that the observed scales crossed the horizon 
during inflation. We call this the pivot scale iV p i vot . All observable quantities are henceforth 
assumed evaluated at this scale. 

The tensor-to-scalar ratio r is defined as the ratio of the tensor power spectrum eq. (3.46) 
to the scalar (curvature) power spectrum eq. (3.47) 

r=jr R ( 3 - 5 °) 

Note that since the tensor power spectrum is conserved on super-horizon scales, the single- 
field result provides an upper bound for the multi-field case. We define the spectral index of 
a power spectrum P x as 

UX = ^k- (3 ' 51) 
Note however that the common definition of the scalar spectral index is related to this 
definition by 

n s = 1 + tir. (3.52) 

We find calculating Trs to be computationally demanding so in practice when calculating 
the spectral index we do the following: 

dlnP x ^ lnP x (N pivot + ^f)- lnP x (N pivot - ^f) _ 
dhik ~ AN 

„ ln P x(Npiyot) ±PX (iVpivot + 1) , 

Px(N pivo t - 1) + iMAWot) ' 

The running of the spectral index a x is a straightforward extension of the spectral 
index, defined as 

a *=dlnT (3 ' 54) 
and hence we approximate it using the same technique, 

~ n x (N pivot + 1/2) - n x (iV pivot - 1/2) . (3.55) 

The minimum number of computations of Trs required to obtain the running in this way is 
three, which is why we take averages when calculating the spectral index. 

We do not consider non-gaussianity in this paper. A methodology for computing it 
within the same formalism has now been given in ref. [17]; however, as for instance shown in 
that paper, it would typically be expected to be small in these types of models. 



4 Findings 

4.1 Two scalar fields 

We performed 5 x 10 6 runs and obtained 75 successful outcomes in terms of sufficient inflation 
achieved without subsequent collapse. Table 3 summarizes some of the mean values accu- 
mulated and shows that our model is in good agreement with observation in all parameters 
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Quantity 


Result 


Observed 


Agreement 


nR 


-0.06 ±0.02 


-0.027 ±0.014 


Y 


OiR 


-0.0003 ± 0.0009 


-0.022 ±0.020 


Y 




0.001 ±0.13 


N/A 


N/A 




-0.02 ±0.22 


N/A 


N/A 


r 


0.05 ± 0.03 


< 0.24 (95% c.l.) 


Y 


Tss 


0.06 ±0.43 


N/A 


N/A 


Trs 


0.8 ±0.9 


N/A 


N/A 



Table 3. Some cosmological parameter constraints. 
Curvature Spectral Index 



Tensor-Scalar Ratio 



-0.08 -0.06 -0.04 

Number of e-folds 



J 



0.30 r 

0.25 

0.20 

0.15 

0.10 

0.05 



0.00 ' f ■ 
0.88 0.90 




0.92 



0.94 0.96 



0.98 



1.00 



1.02 



Figure 2. Histograms of the curvature spectral index, tensor-to-scalar ratio, and e-folds, plus a 
plot of the tensor-to-scalar ratio against the spectral index where the points have been colour-coded 
according to the number of e-folds. The blue line on the n s -r plot is the 95% confidence contour from 
WMAP+BAO+tfo- 



so far tested. Figure 2 shows histograms of the curvature spectral index and tensor-to-scalar 
ratio, a histogram of the number of e-folds and a plot of r against n s , colour-coded with 
the number of e-folds. The last plot also shows the present observational limits from a data 
compilation including WMAP7 results [19]. 

Planck hopes to measure the tensor-to-scalar ratio with an accuracy of a few hundredths, 
hence has discovery potential if it is of order 0.1 or so. For this reason it is interesting to note 
that of the universes not already rejected by the n s -r plot, there seems to be a preference 
for universes with larger r. We also note that all large e-fold universes (green points) lie 
within the 95% confidence limit and only a couple of blue points are rejected, but more data 
are needed to conclude whether or not large e-fold universes are favoured within our model. 
We find this potential trend quite intriguing as it is easy to conceive of measures that give 
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strong weighting to universes with a large number of e-folds. Yet it is generally thought that 
measures exhibiting this kind of favouritism are rejected observationally due to a tendency to 
suffer from what is known as the Q-catastrophe (see ref. [18] for a list of generic problems in 
constructing measures). Our results hint at a tension with this view, since we find universes 
with a large amount of inflation can still agree with observation. Our work does not address 
this problem but it would be interesting to see if this trend is maintained for the very rare 
cases of universes with a huge number of e-foldings, say in the thousands. 

We need to clarify what agreement actually means in this context. Both the observed 
data and the results of our experiment are given in terms of a distribution of values for the 
observables. As it happens, the statistical spread of each is comparable at present, with a 
substantial area of overlap; one can therefore conclude that there is a good chance of a single 
realization from our model giving predictions in accord with the observations. If a model 
data point lies outside the range observed, then it simply indicates that we do not live in that 
universe and such data would not necessarily act against the model. Future observations hope 
to home in on a single value for each observable to ever increasing accuracy, which will clearly 
soon have higher precision than our model predictions. Provided the observationally-favoured 
region remains within the envelope predicted by the model, however, such higher-precision 
measurements will not in themselves be able to argue against the model, at least without 
further model refinements for instance around the choice of measure and anthropic arguments. 
At this point we then must invoke the Copernican principle over this remaining anthropically 
reduced landscape and say we expect to observe values according to the statistics of the model 
distribution. It is in this sense that we say our results are in agreement with observation. 

The histogram of e-folds, Fig. 2, shows that the number of universes drops off rapidly 
with the number of e-folds. We rejected universes with less than 60 e-folds of inflation but 
in the process of doing so we found that sufficient inflation was a very rare process, testing 
millions of universes to find tens of candidates. This is not a new result for "stringy" models. 
For instance ref. [20] investigated one of the most rigorously derived inflationary models 
from string theory, namely brane-anti-brane inflation, and found much the same thing, while 
ref. [12] found the same result for tunnelling landscape models. For our model though this 
is of no concern at all. We are not worried about what proportion of field space allows for 
anthropically suitable conditions, only that there exists some proportion. 

Moving on to the results of the calculations done in the previous section, we look at the 
role played by isocurvature perturbations in modifying the adiabatic perturbations from their 
horizon-crossing value. Figure 3 confirms that the curvature power spectrum always increases, 
and for a non-negligible proportion of the time the fractional change is very large. Looking 
at the histograms for the transfer functions we see that the isocurvature perturbations are 
an important source, since at horizon exit they are of the same order as the curvature power 
spectrum and it is not particularly uncommon for T&s to be significant or even greater than 
one. 

One consequence of the super- horizon fuelling of adiabatic perturbations is that we have 
explicit violation of the consistency relation, as demonstrated in the plot of r against e in 
Fig. 4 (the tensor spectral index in these models is given by nx = — 2e as usual). Substituting 
eqs. (3.46) and (3.47) into eq. (3.50), we have 

16e 

r ~ 1 + ACTrsH/v + T| 
So we see we have a line of points corresponding to r = 16e as in the single scalar field case 
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(4.1) 



Fractional Change in Curvature Power Spectrum 



Fractional Change in Isocurvature Power Spectrum 



8 10 



n , n 



Figure 3. Histograms of the fractional change in the curvature and isocurvature power spectra 
accompanied by histograms for the transfer functions. Arrows indicate Universe 1832942. 



0.14 
0.12 
0.10 
0.08 
0.06 
0.04 
0.02 



Figure 4. Plot of r against e, demonstrating the consistency relation r = — 8rtx = 16e becoming an 
inequality. 



but with many points dropping below the line where Trs and the turn rate have suppressed 
the value of r. Hence the relation becomes an inequality when there are extra degrees of 
freedom, as first described in ref. [21]. 

As an example of the sort of situation that can give rise to large Trs values we include 
Fig. 5. As mentioned earlier, generally the dominant effect is the turn rate and while there 
are a number of types of trajectory that lead to large Trs, in this case we see that it is the 
cumulative effect of a fairly continuously curved trajectory spiralling its way to a minimum. 

As pointed out in ref. [1] (see refs. [19, 22] for primary references), the CMB places strong 
observational constraints against the existence of isocurvature perturbations. As such, any 
multi- field model must either generate no isocurvature perturbations at horizon exit, or must 
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£1 
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Figure 5. Example potential giving rise to Trs = 3.6, and the corresponding evolution in the </>~x 
plane. 

find a way to make them disappear after horizon exit but before horizon entry. Our model 
falls nicely into the second category. Fig. 3 shows that the transfer function T$s consistently 
causes the isocurvature power spectrum to go to zero after horizon exit (except in one case 
which we will discuss in a moment). Referring back to eq. (3.40) and Table 2, we see that 
pure damping occurs when Mmm < and M±± > 0; we find that this scenario is typically 
the case, an example of which is shown in Fig. 6. This result is not particularly surprising 
since one consequence of our choice of measure is that hilltop inflation is a very rare process, 
thus most inflation scenarios occur in the vicinity of a minimum. In such cases, at least 
to some extent the trajectory will be moving towards that minimum. We then have that 
perpendicular curvature will usually be positive and since the minimum must be very close 
to V = 0, the log of the potential will have negative curvature in the direction of motion. 

4.2 The curious case of Universe 1832942 

Figure 3 shows one dramatic exception to the damping of the isocurvature perturbations 
just described. In one case the isocurvature perturbations were actually fuelled rather than 
damped. Fig. 7 shows plots of the potential, its logarithm, the trajectory, and a plot of the 
parallel and perpendicular components of the mass matrix near the end of inflation for this 
case which we call Universe 1832942. Remembering Table 2, we see that for the majority of 
the evolution the isocurvature perturbations are damped but near the end of the evolution 
the trajectory overshoots the minimum causing a sharp reversal in direction as it reaches the 
far side of the valley. As the trajectory goes through this sharp turn at low V the arguments 
above reverse, but we see that the dominant effect in this case seems to be the change in 
curvature as the trajectory traverses the minimum. 

We have severe violation of SRST so our calculations of observables are no longer reli- 
able, but we find the values n s = 0.94 and r = 0.005, with the low r value resulting from a 
significant evolution of the power spectrum after horizon exit of 108%. We see that Universe 
1832942 sits within the 95% confidence contour of Fig. 2, but we remind the reader that this 
confidence limit does not take surviving isocurvature modes into consideration. 




- 18 - 




c— folds 



Figure 6. Example plot of In V, the trajectory, and the relevant mass matrix components, showing 
why Tss tends to be so small. 

4.3 How will observables change with Dl 

As was mentioned at the beginning of this paper, we feel much of the phenomenology of a D 
scalar field landscape is manifest with just two scalar fields. This is true in terms of the types 
of phenomena that exist, but can we make any comments on how the values of observables 
might change with Dl 

The short answer is no. From a statistical standpoint we know that the mean values of 
e and Trs will both increase with D, so one might think that parameters such as n s and r 
could tend to some value as a result of one parameter coming to dominate. This may well 
be true but how so is not obvious. The lack of any strong correlation in Fig. 8 for n s with 
respect to e or Trs leads us to think that predicting how n s will change with D is non-trivial. 
That said, there is an important distinction between how e and Trs will change with D. 
Trs is unbounded, in that as you increase D you just get an increasing sum of Tr5 4 terms; 
this has both a cumulative effect and increases the statistical chance of encountering a large 
TrSz ■ The value of e in contrast is capped, if for no other reason because inflation will end if 
e is too large, but it may not ever even increase that much. 

There are a number of examples in the literature where a large number of scalar fields 
can actually lead to a decrease in e, one example being assisted inflation [23]. So for cases 
where the absolute values of e and Trs are important we can hazard a guess and say that 
for sufficiently large D eventually the dominant effect will be Trs- This is the case for the 
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Figure 7. Plots of the potential, log of the potential, trajectory and relevant components of the mass 
matrix (M\\\\ is blue, M±± is red) for the one case we found, Universe 1832942, where the isocurvature 
perturbations were fuelled rather than going to zero. The arrow indicates where the trajectory turns 
back on itself. 



1.5 2.0 



Figure 8. Correlations of n s with e and Trs. 



tensor-to-scalar ratio and remembering Fig. 4 and eq. (4.1) we would thus expect increasingly 
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severe violation of the consistency relation, with a suppression of r as we move to higher D. 
We will investigate this as part of a future study of the D-dimensional case. 

5 Discussion 

We have investigated the properties of inflationary trajectories in a toy-model landscape 
with two scalar fields, in which we fully track the effects of isocurvature perturbations. We 
have focussed on the well-motivated case where the variations in the potential correspond to 
masses of the order of Mp\. Trajectories with sufficient inflation are rare, with one successful 
run per roughly 10 5 randomly generated potentials, but the successful runs are typically in 
good accord with observational constraints on the perturbations, mostly lying within the 
region of the n s -r plane delineated by the WMAP+BAO+ffo 95% confidence contour. 

We find that isocurvature perturbations naturally go to zero after horizon exit. This is 
a direct result of the typical geometry of the landscape in the vicinity of an inflationary tra- 
jectory. For isocurvature perturbations to not go to zero one requires an unstable trajectory, 
which happened only in one of our successful realizations. Nevertheless, in many cases the 
isocurvature perturbations have a lasting consequence as we find significant fueling of the 
adiabatic perturbations from the entropy perturbations to be a fairly common occurrence. 
Because of this, we find that the tensor-to-scalar amplitude often lies below the level that 
would be predicted by the single-field consistency equation. One might speculate that this 
effect is likely to become more prominent as the number of scalar fields is increased, both 
from a statistical standpoint and as a purely cumulative effect, since such evolution is closely 
linked with turns in the trajectory. We plan to investigate this dependence in future work. 

While our work significantly extends previous comparable analyses in the literature, 
such as ref. [4], it nevertheless remains rudimentary and many further steps are needed 
if it is to become more realistic, beyond the issue of the number of scalar fields. Most 
pressing is a proper treatment of the way the measure problem affects our work, particularly 
in combination with anthropic selection. At the moment we have largely ignored this and 
obtained results which are based on a uniform distribution of initial conditions, rather than 
as drawn from random late-time observers. 
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